A causal look into the quantum Talbot effect 
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A well-known phenomenon in both optics and quantum mechanics is the so-called Talbot effect. 
This near field interference effect arises when infinitely periodic diffracting structures or gratings 
are illuminated by highly coherent light or particle beams. Typical diffraction patterns known 
as quantum carpets are then observed. Here the authors provide an insightful picture of this 
nonlocal phenomenon as well as its classical limit in terms of Bohmian mechanics, also showing the 
causal reasons and conditions that explain its appearance. As an illustration, theoretical results 
obtained from diffraction of thermal He atoms by both A^-slit arrays and weak corrugated surfaces 
are analyzed and discussed. Moreover, the authors also explain in terms of what they call the Talbot- 
Beeby effect how realistic interaction potentials induce shifts and distortions in the corresponding 
quantum carpets. 

PACS numbers: 



I. INTRODUCTION 

Particle diffraction by different types of devices (slits, 
gratings, or surfaces) has become a standard technique 
to test the validity of quantum mechanics. This is con- 
firmed by a large amount of experiments, ranging from 
tiny objects (e.g., electrons, neutrons, single atoms, or 
small clusters) to more complex, mesoscopic-size systems 
(e.g., fuUerenes, large biomolecues, or Bose-Einstein con- 
densates). Experiments with large, structured particles 
induce to think that, in principle, there are no size- 
dependent limits to observe particle diffraction, except 
those associated with the preservation of the system co- 
herence. This is a very important issue, since the system 
coherence may disappear relatively fast as a consequence 
of the many eventual degrees of freedom involved in the 
process, either related to the surrounding environment 
(external) or associated with the own structure of the 
particle (internal). 

In 1836, when characterizing optical gratings, Talbot^ 
observed a repetition of alternate color bands of com- 
plementary colors (red-green and blue-yellow) at certain 
distances from the grating. About 50 years later, in 
1881, Rayleigh^ proved that this phenomenon is a con- 
sequence of the diffraction of highly spatially coherent 
(plane) waves by gratings; the color band structure ob- 
served is a manifestation of the periodicity and the shape 
of the grating. The alternation of complementary color 
bands occurs at integer multiples of Zt = d^/A, the Tal- 
bot distance^ where d and A are the grating period and 
the wavelength of the incident plane wave, respectively; 
bands with equal color patterns thus repeat at integer 
multiples of 2zt- The Talbot effect has important tech- 
nological applications in optics, such as image processing 
and testing or production of optical elements. Similarly, 
its quantum-mechanical counterpart is relevant in elec- 
tron optics, where it has many applications within elec- 
tron microscopy. Moreover, this effect has been observed 
experimentally with heavy particles, such as Na atoms- 



or Bose-Einstein condensates.— 

Like in optics, A^-slit arrays can also be considered typ- 
ical examples of quantum gratings. When such arrays are 
illuminated by continuous, coherent wave fronts, a con- 
tinuous quantum flow is observed behind the slits. This 
flow displays a typical pattern called quantum carped 
with periodicity d along the direction parallel to the plane 
containing the slits {x) and 2zt along the propagation 
direction (z). Full recurrences along z are found at in- 
teger multiples of zt] these recurrences are the direct 
analog of the color bands observed when working with 
optical gratings. A recurrence that coincides with the 
initial state describing the (diffracted) system is called a 
revival of such a state and appears at integer multiples 
of zt\ for even integers, the state looks exactly the same 
as the initial one, while for odd integers it is shifted half 
a period (c?/2) with respect to the latter. Recurrences at 
fractions of zt consist of superposed images of the initial 
state with itself. Indeed, if the boundary conditions of 
the slits are "sharp" (the window function is not differen- 
tiable at the borders of the slit), one can observe fractal 
structures at irrational fractions of zt^ which give rise to 
fractal carpets ^^^^2^^ 

The Talbot carpet associated with the diffraction of 
a monochromatic beam of He atoms [E^ = 21 meV, 
A = 2'nhl^J2mE:, = 0.991 A) by a grating consisting 
of 50 Gaussian slits with period d — 3.6 A is displayed 
in Fig. [TJa). Here, the x variable is normalized to the 
grating period, d, and the z variable to twice the Talbot 
distance 2zt, with zt = 13.08 A (throughout this work 
we will follow the same convention). The carpet appears 
when we record pt{x,z) = |5'j(a;,2;)p along x at differ- 
ent times, monitored in terms of z (the wave propagates 
along the z direction at a constant speed Vz — 27r/i/mA; 
see below for details). That is, at a given time ipiot we 
have plotted the "slice" Ptpiot (2;, Zpiot) of the probabil- 
ity density, where x runs over all possible values along 
a certain range (in the figure, from —6 to 6) and the z 
value is fixed: Zpiot — ^(ipiot) = t^z^piot- Some of these 



2 



, .{i] 



1^ 




X 




z = 



= V4 



z = z /2 



z = 3z/4 



FIG. 1: (Color online) (a) Quantum carpet for He diffracted 
by a 50-slit array. The color scale, from blue to red, indicates 
increasing values of the probability density, (b) Snapshots of 
the probability density from z — to z = zt- In both panels, 
the X distance is scaled in units of the grating period d and z 
in units of twice the Talbot distance {2zt = 2d^ /X). 



of the grating size. At such distances the carpet gradually 
fades, and the typical Fraunhofer patterns that character- 
ize most particle diffraction experiments start becoming 
apparent. 

Although Talbot patterns are well known in the liter- 
ature, as far as we know there is no detailed explanation 
for them in causal terms. This work is aimed to pro- 
vide a full causal interpretation for this phenomenon as 
well as for its relationship with the Fraunhofer diffraction 
and the classical limit. For this purpose, we have chosen 
Bohmian mechanics, where the standard wave picture is 
replaced by trajectories in configuration space. The use 
of such trajectories is very important in order to shed 
some light into real experiments, where a picture in terms 
of the motion of individual particles is always highly de- 
sirable (specially if particles follow the quantum flow, un- 
like other classical representations). This is an important 
issue, for instance, in surface science experiments, where 
the illumination of a surface for its study and characteri- 
zation is limited in extension. Note that, when assuming 
that no imperfections are present, surfaces can be con- 
sidered as gratings analogous to iV-slit arrays because of 
their periodicity. Nonetheless, the short-range attractive 
interaction undergone by the diffracted particles near the 
surface leads to a slight shift and distortion of the Tal- 
bot period. We have called this effect the Talbot-Beeby 
effect since the variation of the diffracted beam wave- 
length due to presence of attractive wells is known as the 
Beeby correctior}^ in surface scattering. 

This work is organized as follows. In order to be self- 
contained, a rigorous mathematical approach in terms of 
standard quantum mechanics will be presented in Sec. [Til 
analyzing both the appearance of the Talbot effect in 
gratings of Gaussian slits and the transition to the Fraun- 
hofer regime. In Sec. IIIII we describe the application of 
Bohmian mechanics to the problems discussed in the pre- 
vious section. A discussion of the Talbot and Talbot- 
Beeby effects in diffraction of He atoms by both an A^- 
slit array and the Cu(llO) surface is given in Sees. IIVI 
and |Vl respectively. In the light of these results, the 
meaning of the classical limit and the quantum-classical 
correspondence is discussed in Sec. lVIl Finally, the main 
conclusions from this work are summarized in Sec. IVIII 



II. WAVE APPROACH TO THE TALBOT 
EFFECT 



slices are shown individually in Fig. [TJb), from z = to 
z = Zt- As mentioned above, revivals of the initial state 
occur at zt (shifted half a period) , while the other slices 
manifest vague reminiscences of ^0(2^,0). After reaching 
z = Zt, the diffracted beam will evolve towards 2zt fol- 
lowing the inverse sequence to that shown in Fig. [TJ^b) 
[see Fig. [UJa)], then an exact revival of the initial proba- 
bility distribution will be observed. This process repeats 
endlessly unless the grating has a finite size. With real 
gratings its extension is limited to distances of the order 



Rigorous analytical studies of the Talbot effect from an 
optical viewpoint can be found, for instance, in Ref. [ill . 
Here we provide an alternative quantum-mechanical 
derivation highlighting the physical aspects underlying 
this phenomenon, in particular, (i) the role of the super- 
position principle and (ii) the analogy/equivalence be- 
tween Talbot patterns and the carpets observed in multi- 
mode cavities (e.g., waveguides). For simplicity, we focus 
on gratings constituted by Gaussian slitSfi^ii^ i.e., char- 
acterized by Gaussian transmission functions, though our 
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analysis can be straightforwardly generalized to any kind 
of periodic grating. Gaussian transmissions can be ob- 
served, for instance, when studying the diffraction of a 
monochromatic beam (with wavelength A) by a soft, re- 
pulsive (exponential) potential barrier with an infinity of 
identical holes (slits) 



1. Single Gaussian slit 

To understand the physics associated with the time 
evolution of the diffracted wave function by an infinite 
Gaussian grating, it is important to consider first the 
dynamics of a single diffracted Gaussian function. As- 
suming the initial time (t = 0) as the instant when the 
incident wave has just passed through the slit with per- 
pendicular velocity = 0), the initial state can be ex- 
pressed as 



-x'^ /4al-z^ /4al+ip^z/h 



(1) 



where A{0) = is the norm and E^nd 

the initial widths in the x and z directions, respectively. 
Within this description, note that Ox is related to the 
transmission function of the slit, giving an idea of the 
aperture dimension (i.e., the slit width); on the other 
hand, cr^ is an indicator of the monochromaticity of the 
incident beam (the limit ct^ — s- cx) represents the case of 
a fully monochromatic incident wave). 
The time propagation of $o is given by 

$t(a:, z) = A(t) e-^V45x.t<T.-(2-20V4*.,t<T. 
with A{t) = {2'Kax,tcrz.t)~^^'^, zt = v^t, ^pl/2m, and 



i = x,z, (3) 



where 




X, z 



(4) 



gives the instantaneous width along each direction. From 
Eq. @ a time scale tg — 2maf /h separating two regimes 
with different "spreading rates" can be defined. If i ^ 
tg, the Gaussian widths almost remain the same {at^t ~ 
Ci). On the other hand, if i ^ tg, the Gaussian widths 
undergo a linear increase with time ((Ti,t ~ ht/2mai). In 
this way, choosing (Jx — d/8 and az = d, $t will not 
almost display any increase in its size along z, but only 
along X. Accordingly, Eq. ([2]) can be reexpressed as 



^t{x,z) = A'{t)e^ 



X e 



ipz{z — zt)/h+iEzt/h 



where A'(t) = {2'Kdx,t'yz)~^^^ ■ Note that this will en- 
sure the overlapping of Gaussian functions (along the x 
direction) in the case of the grating, and therefore the 
appearance of interference. 

For hard slits, the interaction potential is zero every- 
where except along the grating (2 = 0), where it is infin- 
ity. Thus, there is no coupling between both coordinates, 
X and 2, and the time evolution of the wave function 
along each direction can be studied separately. Since the 
interference takes place along x, we will only focus on the 
part of the wave function depending on this variable. 



0t(a;)=^,(i)e-^'/4a.>,.. 



(6) 



where Ax{t) = (27r(T^ Regarding the component 

of the wave function along z, it is only important to know 
that it propagates at a constant speed Vz, i.e., its cen- 
troid evolves according to z{t) = v^t. This allows to sim- 
plify our analysis, reducing it to a one-dimensional time- 
dependent problem, where the z coordinate contains the 
same information as t. Thus, from now on the subscript 
X will be dropped from the magnitudes related to the 
spreading along the x direction (cr and A instead of Ux 
and Ax^ respectively), and from the momentum associ- 
ated with the X coordinate (p instead oi Px)- 



2. Infinite periodic slit gratings 

According to Bloch's theorem;^ the problem of finding 
the diffracted wave function associated with a given infi- 
nite periodic potential reduces to determining the wave 
function associated with a single unit cell of such a po- 
tential; the full wave function is just a repetition of the 
latter, which satisfies the following Born-von Karman 
boundary conditions: 



iJtix + d) = ijjtix). 



(7) 



Thus, we assume that the total number of Gaussian 
slits amounts to = 2K + 1, the slits being centered 



q*^"* = kd (with k 



0, ±1, ±2, 



±^0 and 



0. Our problem then reduces to only consider 
one of the Gaussian wave packets (pt, constituting the 
full wave function. We can choose, for instance, (pt to 
be the wave packet with fc = 0, which is confined within 
the space region x = ztd/2 and whose time evolution 
is described by Eq. In direct analogy to surface 

science, from now on we will call the space region of 
length d enclosing each (initial) Gaussian wave packet 
[i.e., {k - l)d/2 <x< kd/2\ a unit cell. 

As is well known, any unbound wave function can be 
represented as a superposition of plane waves. 



with 



, (5) 



4>t{x) 



a{p) 



^ipx I h—iiut 



dp, 



<i)o{x)e-'P''/^dx. 



(8) 



(9) 
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Thus, Eq. dni) can be expressed as 

(8^^2)1/4 



^-a^p^ /h^+ipx/h-iujt I'-i^Q-j 



Because of the periodic boundary condition ([7|), not all 
momenta are allowed; we then pass from a continuous 
basis of momenta to a discrete one, and Eq. pO|) becomes 



2\ 1/4 




<t>t{x) 



(11) 

where p„ — 2Trhn/d and ojn = 2Tr'^ hn'^ / mcP . 

The initial (full) wave function describing the system 
(along the x direction) is assumed to be a coherent, 
nonoverlapping superposition of identical Gaussian wave 
functions that propagate along the z direction with con- 
stant velocity. 



K 



ipo{x) oc lim A{0) < 

k=-K 



-{x-kdf/4a^ 



(12) 



The corresponding time-evolved wave function is 
straightforwardly obtained by replacing the Gaussian 
wave packets in Eq. (|12p by their time-dependent coun- 
terparts given by Eq. pTjl . This leads to 



K—*oo 



(13) 



after rearranging terms. As one would expect, Eq. ([Tc 
indicates that all the information regarding the infinite 
grating is contained within a single unit cell, with the 
factor 2K + 1 arising from the total number of unit cells 
considered. To avoid the divergence introduced by this 
factor, we reexpress Eq. as 



lim 

K^oc 



A{t) 



K 



-{x~kd) /4ata 



2K 



1 ^ 

k=-K 



(14) 



The smallest time elapsed necessary to observe a recur- 
rence of the wave function is determined by the condition 



(15) 



Using Eq. pT|) . it is straightforward to see that this con- 
dition is satisfied whenever e"*'^"'^'' — e~^'^' for all n, i.e., 
when Tr = 2Tr/LUi = uKp/Trh. At this time, the z dis- 
tance between two consecutive recurrences of (j)t will be 
Zr = VzTr = 2cf / X, i.e., twice the Talbot distance (2zt)- 
As shown in Fig. [Ijb) , there are also recurrences with pe- 
riodicity Zr — cP/X — zt- These recurrences result from 
considering the symmetry 



A+Tj2{x + d/2) = 'ipt{x), 



(16) 



which appears when taking the terms e^Vnxlh luj^r^ ^ 
whole. 




FIG. 2: (Color online) Schematic picture of the formation of 
revivals in (a) gratings of Gaussian slits and (b) multimode 
cavities. Arrows indicate the direction of the motion of the 
"interfering" disks (see text for details) and straight lines rep- 
resent the corresponding paths. In part (a), the numbers label 
different slits, and the blue line limits the Talbot region. 



The observation of revivals can be explained by means 
of a simple schematic picture, which is displayed in 
Fig. [IKa). To understand this picture we consider that 
each wave packet can be represented by two "interfering" 
hard disks that propagate along the z direction with op- 
posite X velocities (the x motion of the center of mass 
is null, as expected from the propagation of a real wave 
packet, which only advances along z). These disks thus 
indicate the spreading of the wave packet. Following the 
motion of the disks associated with A; = 0, we observe 
that they meet at a; = ±d/2 and z — zt with those 
arriving from the neighboring slits (fc = ±1). Since all 
disks are identical, there will be full constructive "inter- 
ference" , this leading to the first revival. The next revival 
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will occur at z = 2zt and is caused by the interference of 
the disks coming from k = ±1, but not from fc = 0; the 
disks leaving the slit with fc = will contribute to the 
revivals observed aX x — ±d and z — 2zt. 

From Fig. [UJa) we can also infer an effective size for 
Talbot structures when the grating periodicity is limited. 
In the picture we have represented five slits; thus, after 
the ingoing-moving disks corresponding to A: = ±2 have 
interfered, no Talbot revivals will be observed. From the 
point of view of Gaussian wave packets, this will happen 
after the size of those corresponding to the outermost shts 
is large enough as to make them to interfere, i.e., when 
2<Tt ~ 2Kd. Assuming that at that stage the width of 
the Gaussian functions increases linearly with time, i.e., 
at ~ ht/2ma, the maximum time at which a (partial) 
revival can be observed is imax ~ 2Kdma/h. Since the 
propagation along z is also linear with time (at the speed 
Vz), the z distance where such a revival will be observable 
is 



Vztn 



2zT 



7r(iV- 1) 



(17) 



(we have particularized this expression to our case, where 
a = d/8). Examples of the validity of this relation will 
be seen in Sec. IIVI Beyond Zmax a transition towards 
the Fraunhofer diffraction regime, with its characteristic 
fringe patterns, starts to be observed. This transition 
takes place beyond the boundaries of the Talbot region, 
enclosed in Fig. HJa) by a blue triangle with height Zmax 
and basis {N — l)d. 

As soon as the wave function (or part of it) leaves the 
Talbot region, Bloch's theorem (and therefore the Born- 
von Karman boundary conditions) is no longer applica- 
ble. To describe the Fraunhofer diffraction, one has to 
start from Eq. (|14|) . but having in mind that the grating 
has a finite size. Thus, at a relatively long z distance from 
the grating, its dimensions are negligible when compared 
with the X distances involved in the diffraction process 
[i.e., x:s> {N — l)d], and therefore 



-{x — kd)'/4at'^ 



-ikKdx / z 



(18) 



where k, = 27r/A. Substituting the right hand side of this 
expression into Eq. (fT4)) we obtain 



V't(a;) 



2K+1 



Z ^iKX IZ 



/2z 



K 



-ikndx I z 



(19) 

Now, using the so-called paraxial approximation from 
opticsiS. and considering xjz = tan 6 « sin 6* {9 is the 
observation angle), the probability density (diffraction 
pattern) along x can be expressed as 



gt{x) = \A{t)\' e- 



siii{NK,dsme/2) 
TV sin {nd sin 0/2) 



(20) 



In this expression, the term between square brackets is 
the structure factor, which accounts for the interference 



among the different diffracted (Gaussian) beams. On the 
other hand, the normalized exponential is the form fac- 
tor, which is related to the diffraction by a single unit 
cell. Because of the information provided by these fac- 
tors, they are very useful to characterize optical gridsi^ 
as well as periodic surfaces?^ 

From the structure factor we note that the Fraunhofer 
fringes appear in accordance to the diffraction (quanti- 
zation) condition 



sin 9 



A 



£ = 0,±1,±2,.. 



(21) 



where £ is the diffraction order. On the other hand, there 
will be a series of minima whenever 



sin0 



f_ A 



\e' -l\ = l,2,...,N-l. (22) 



As seen, there are iV— 1 minima between two consecu- 
tive principal maxima and, consequently, iV— 2 secondary 
maxima; the height of the latter is smaller than that of 
the principal maxima and decreasea^^ very fast with N. 

Finally, a comment regarding the local spreading of 
each Gaussian wave packet and the boundary condition 
Eq. ([7]) is worth stressing. According to Eq. ([6]), after a 
certain time the size of every Gaussian in (|14p at t = 
will be such that it will extend beyond the boundaries of 
its corresponding unit cell. That is, the spreading of (l)t, 
which is a local property, gives rise to the appearance of a 
nonlocal behavior, where each part of the resulting wave 
function, ipt, is strongly influenced by the presence of the 
remaining ones. Note, indeed, that the Talbot effect is 
precisely a nonlocal effect: it emerges as a consequence of 
the overlapping of many identical wave packets. As will 
be seen below, Bohmian mechanics provides a natural 
picture to this problem in terms of the regimes of motion 
associated with two well-defined dynamic equilibria. 



3. Multimode cavities 

The Talbot effect is closely related to multimode 
interferencefi^iiiii^ the interference process that takes 
place when an infinity of modes of a cavity are super- 
posed. Waveguides are a typical example of multimode 
cavity where the wave motion is constrained along one 
direction (e.g., x) and unbound along the other (e.g., z). 
Assuming no coupling between the motion along each 
direction, and that the cavity is a square box along x, 
centered at a; = and with length d, the time evolution 
of a Gaussian wave packet can be expressed at any time 
in terms of the modes corresponding to this cavity as 



= V 3 



d\ 

OQ 

e-'^'p"^/^"-''^-*/^ cos{pnx/h), n = 0, 1,2, . 



n=0 



(23) 
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with pn = {2n + l)7rh/d and £"„ = p\/2m. It is easy to 
show^^ that recurrences in the probabihty density arising 
from Eq. ([23]) occur at integer multiples of the period 
associated with the smallest frequency. In our case, this 
frequency is a;i,o — {Ei — Eo)/h = ATr^h/mcP, and the 
period associated is 



27r 
Win 



md 
27rft' 



(24) 



We have to mention that the same periodicity can also 
be observed in the wave function, except for a constant 
phase factor [(pt+Tr.{x) = e''^0t(x)]. 

Since the wave function evolves at a constat speed 
along z, the recurrences will also give rise to the for- 
mation of quantum carpets; some examples of these car- 
pets can be found, for instance, in Ref. [H. These re- 
currences appear at integer multiples of the distance 
Zr — VzTr = d^/A, which corresponds to the Talbot dis- 
tance obtained in the case of a periodic grating. This is 
the analogous case to the revivals observed at z = 2zj' 
in grating systems. To understand why the same type of 
recurrences cannot be observed in both periodic gratings 
and cavities, it is very instructive to look at the scheme 
represented in Fig. [IJb). This picture is equivalent to 
that shown in Fig. [5)^ a), but with the difference that the 
only possibility for the disks to cross is when they meet 
again; no interference at a; = ±d/2 is possible. In other 
words, while the quantization leading to Eq. pip arises 
from a "matching" condition at the borders of neighbor- 
ing unit cells (thus allowing interference), in Eq. it 
comes from having impenetrable boundaries. 



III. BOHMIAN APPROACH TO THE TALBOT 
EFFECT 

The fundamental equations of Bohmian mechanics are 
commonly derived by expressing the system wave func- 
tion in polar form^i^iiS^ 



«'t(a;,z) =py^(a;,z)e'^*(^'^)/'^. 



(25) 



Here, pt is the probability density and St is the (real- 
valued) phase. Substituting Eq. (|25p into the time- 
dependent Schrodinger equation and making a bit of al- 
gebra, one reaches two (real-valued) coupled differential 
equations 



dSt (WSt)^ 
dt 2m 



0, 



V + Qt = Q. 



(26) 
(27) 



Equation (|26p is a continuity equation that ensures the 
conservation of the quantum particle flux. On the other 
hand, Eq. (j27p . more interesting from a dynamical view- 
point, is a quantum Hamilton- Jacobi equation governing 
the particle motion under the action of a total effective 



potential V^^ = V + Qt- The last term in the left hand 
side of this equation is the so-called quantum potential 



2 1/2 



Qt 



Pt 



(28) 



This context-dependent, nonlocal potential determines 
together with V the total force acting over the system. 

In the classical Hamilton- Jacobi theory {Qt = 0), St 
represents the action of the system at a time t, and the 
trajectories describing the evolution of the system corre- 
spond to the paths perpendicular to the constant-action 
surfaces at each time. Similarly, since the Schrodinger 
equation can be rewritten in terms of the Hamilton- 
Jacobi equation ([?f|) . St can be interpreted as a quan- 
tum action satisfying similar mathematical requirements 
as its classical homologous. In Bohmian mechanics the 
classical-like concept of trajectory thus emerges in a nat- 
ural manner; particle trajectories are defined as the so- 
lutions of the equation of motion 



VSt h r t 1 

i = — Im ^-r V*t 

mm 



(29) 



Moreover, since in Bohmian mechanics the system con- 
sists of a wave and a particle, it is not necessary to 
specify the initial momentum for each particle (as hap- 
pens in classical mechanics) , but only its initial position, 
To = (s^Oj^o); a-s well as the initial wave function ^'q. 
The initial momentum field is predetermined by ^'o via 
Eq. ([29]) . and the statistical predictions of the standard 
quantum mechanics are reproduced by considering an en- 
semble of (noninteracting21) particles distributed accord- 
ing to the initial probability density po = l^oP- 

Equation (^5)1 is well defined provided the wave func- 
tion is continuous and differentiable. This is not the case, 
however, for quantum fractals,'^'^-*'?' where the Bohmian 
mechanics based on Eq. (|29p is unable to offer a trajec- 
tory picture for the corresponding wave functions. This 
apparent incompleteness can be nevertheless "bridged" 
by taking into account the decomposition of the quan- 
tum fractal as a sum of (differentiable) eigenvectors of the 
Hamiltonian, and then redefining Eq. (|29p conveniently. 
This is seen in detail in Ref. [H, where a generalization 
of the standard formulation of Bohmian mechanics to in- 
clude quantum fractals is given. 

Prior to any calculation, some physical insight into slit 
diffraction problems is possible by studying the proper- 
ties of the velocity field. Indeed, the fact that Bloch's 
theorem (together with the Born-von Karman bound- 
ary conditions) holds simplifies this study from both a 
conceptual and a computational perspective. Concep- 
tually, because the study of the full system reduces to 
only understanding the dynamics within a single unit 
cell. In this sense, the analysis is similar to that of hav- 
ing a multimode cavity, as said above. This implies a 
computational advantage: it allows to perform calcula- 
tions taking into account only what happens inside a 
single unit cell and periodic boundary conditions [i.e., 
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ipt{x—d/2) = ipt{x+d/2)], thus reducing the computation 
efforts. Of course, this simphfication is only possible if 
we are under the assumption of infinite periodic gratings 
or, at least, we are working well inside the Talbot area, 
delimited by the blue triangle shown in Fig. [2ja). Other- 



wise, this advantageous framework is no longer valid and 
the full system wave function has to be considered. 

According to the previous statements, substituting 
Eq. pT|) into ((29l) for the x coordinate yields 



1 , K e-'(p'+P?)/''' cos[(p, - pj)x/h - (iu, - uj,)t] 



m 



C0s[(pi — Pj)x/h — {uJi — LUj)t] 



i,j-±l,±2,. 



(30) 



From this equation we can extract relevant information 
about the physical properties of the quantum trajectories 
(or, equivalently, their topology). Note that the velocity 
field ((30)) satisfies exactly the same symmetry conditions 
expressed by Eqs. and The particle motion is 

thus oscillatory, with the recurrences displayed by pt oc- 
curring in those space regions where trajectories accumu- 
late. However, unlike the disks depicted in Fig.[2l^a), the 
trajectories associated with each single unit cell will re- 
main inside it at any time and will never cross those com- 
ing from other unit cells: the particle motion is bound. 
This fact is a manifestation of the noncrossing property of 
Bohmian mechanics: trajectories can never pass through 
the same point on configuration space at the same time 
due to the single valuedness of the momentum field. This 
explains why the physics along each unit cell (or, by ex- 
tension, the physics associated with infinite gratings) is 
so close to that observed in multimode cavities. The 
boundary periodic conditions give rise to the presence of 
nonphysical impenetrable walls at x — ±c?/2, where the 
velocity field ([50]) vanishes and therefore the correspond- 
ing quantum trajectories will be just straight lines. Due 
to the noncrossing property, these trajectories will act as 
(impenetrable) boundaries for particles coming from dif- 
ferent neighboring slits. This is the Bohmian causal ex- 
planation for the Born-von Karman periodic boundary 
conditions. The same effect is also found for the trajec- 
tories starting at the center of the slits (a; = 0) , since the 
velocity field is also zero along this symmetry line. These 
trajectories evolve along two directions which are spec- 
ular to one another (with respect to a; — 0). Note that 
this goes beyond the classical-like picture provided by 
the disks associated with a multimode cavity in Fig.[2Jb) 
in the sense that it adds a constraint of different nature 
to the types of motion that one can observe in quantum 
mechanics. Moreover, this behavior is regardless of the 
motion along the z direction (provided there is no cou- 
pling between both directions), where the corresponding 
equation of motion renders 

z{t)^zo + v,t. (31) 

That is, all particles display the same uniform rectilin- 
ear motion, as happens in classical mechanics, since the 
swarm of particles is basically guided by a plane wave. 



If a size-limited grating is considered, after some time 
the trajectories will be out of the Talbot region, and 
therefore the description given above will no longer be 
applicable. Far beyond the grating, one can appeal to 
the Fraunhofer approximation in order to gain some in- 
sight on the topology of the trajectories. In such a case, 
introducing Eq. into (|^ leads to 



Assuming that the probability density is only significant 
along the quantized values oi x/z w sin 9, where one 
observes the Fraunhofer principal maxima (we neglect 
the presence of secondary maxima), Eq. ([5^ becomes 

x^iv,^. (33) 

This assumption is equivalent to considering that the 
diffracted wave function consists of different independent 
plane waves, each one characterized by a quantized mo- 
mentum p£ = £{2TTh/d). This is a very important result 
that can be understood in terms of two different equi- 
librium regimes. The first equilibrium regime occurs in 
the Talbot/Fresnel region and is characterized by what 
we could call an equilibrium of momenta. That is, in this 
region the possible momenta satisfying a certain quanti- 
zation condition are selected. This selection depends on 
the features defining the unit cell of the grating (or, in 
other words, the interaction potential between the grat- 
ing and the diffracted particles). The second equilib- 
rium regime, which we could call equilibrium of configu- 
ration, happens far beyond the grating, once the particle 
distribution remains with the same shape (regardless of 
spreading effects) . The transition from the momentum to 
the configuration equilibrium is a direct consequence of 
the redistribution of momenta among the different parti- 
cles contributing to pt . These momenta make the swarm 
of particles to evolve in such a manner that, at a certain 
distance from the grating, it will separate in different 
beams, each moving with a different (average) momen- 
tum Pi — £{27Th/d). The equation of motion for the par- 
ticles will then be 

27rM 

xit) ^xo + — J- 1, (34) 
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where £ = denotes the classical direction of motion. 
This contrasts with the remaining beams, which undergo 
a classical-like motion though they do not follow the 
real classical trajectory; the residual term involving ft 
is a clear indicator of the nonlocal behavior of quantum 
mechanics.— Nevertheless, as m increases, this deviation 
from the real classical motion gets smaller and smaller 
but never zero (as will be shown below). 

The effects derived from the noncrossing property in 
periodic grating systems are equivalent to considering the 
presence of some effective impenetrable barriers. In the 
same way, the temporary accumulation of quantum tra- 
jectories in certain space regions could also be thought as 
the effect of an effective quantum pressure. For instance, 
the particles moving along an impenetrable (effective) 
barrier will keep their motion until those arriving from 
the central part of the initial wave packet will move to- 
wards a; = 0, thus decreasing the "pressure" exerted on 
the former (the same effect has been discussed in multi- 
mode cavitiesi^ or soft double slit&i^). The decrease of 
the quantum pressure is responsible for the appearance 
of Fraunhofer fringes; beyond the Talbot region, the par- 
ticles feel a smaller pressure that allows them to get into 
different channels. After some time, particles reach the 
configuration equilibrium regime, where they move along 
well-defined channels, the well-known Bragg diffraction 
channels. 



IV. SLIT ARRAYS 




First we consider the dynamics due to slit arrays, 
where He atoms are the diffracted particles and the slits 
are described by Gaussian (transmission) functionsj ^^i^^ 
In Fig. [2Ia) we can observe that, in accordance to the 
previous discussion in Sec. lIIIi trajectories follow the flow 
characterizing pt in Fig.[T](a). Note that here the trajec- 
tory carpet is not a one-piece structure, as happened with 
the Talbot carpet generated by pt in Fig. [T](a) , but it con- 
sists of many unit substructures (as many as slits have 
been considered). These structures are the causal effect 
of having indistinguishable unit cells. Moreover, trajecto- 
ries exiting from one of the slits never cross those leaving 
the other slits, this being due to the noncrossing property 
described above. This can be better seen by looking at 
the enlargement presented in Fig. [3]Jb) . In this plot we 
observe that initially the trajectories leave the slit in a 
diffusive manner towards each border of the slit. Since all 
the slits are identical (as well as their transmission func- 
tion), these trajectories will feel in a short term the pres- 
ence of bunches of trajectories coming from the neighbor- 
ing slits. Then, the trajectories start bending until they 
move perpendicular to the slits for a while, pushed by the 
neighboring trajectories. This is a clear manifestation of 
the quantum pressure: the pressure exerted by bunches 
of trajectories moving in opposite directions gives rise to 
an effect similar to that of having an impenetrable (infi- 
nite) potential barrier. But the quantum pressure is also 



FIG. 3: (Color online) (a) Quantum trajectories associated 
with Fig. [ija). Trajectories in different color indicate that 
though according to Fig. [TJa) the Talbot carpet seems to be 
a one-piece pattern, it is indeed constituted by many sin- 
gle bunches of trajectories arising from each slit and without 
crossing with those exiting from neighboring slits, (b) En- 
largement of part (a) in the region close to the slits in order 
to show the correspondence with Fig. \Vlh) . In both panels, 
the X distance is scaled in units of the grating period d, and 
z in units of twice the Talbot distance {2zt)- 



felt from the action of the trajectories coming from the 
same slit: those started with initial conditions closer to 
the center of the slit will push the other to move parallel; 
only when they start moving back again towards the cen- 
tral axis of the slit, the quantum pressure will decrease 
enough as to allow the outer trajectories to move again 
towards the original position of the slit. In this way, at 
z = 2zt we will recover again the initial pattern of a sum 
of Gaussians. It is interesting to stress that the max- 
ima at z = zt, shown in Fig. [H^b), have not the same 
structure as those at z = or z = 2zt; the trajectories 
contributing to the maxima at z — zt start from two 
different (neighboring) slits, while a,t z = 2zt only the 
trajectories started from the same slit will contribute to 
the corresponding maximum. 

Provided the grating extends to infinity, the unit struc- 
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FIG. 4: (Color online) Left: Appearance of the Talbot carpet 
within a certain space region as the number of slits increases: 
(a) iV = 1, (b) iV = 2, (c) = 3, (d) iV = 10, and (e) 
N = 50. Right: Quantum trajectories corresponding to the 
cases shown in the left panels. In all panels, the x distance is 
scaled in units of the grating period d, and z in units of twice 
the Talbot distance (2zt). 



tures of the pattern seen in Fig. ^a,) repeat indefinitely. 
However, standard slit arrays have a limited size as well 
as the the incident beam. Therefore, it is clear that 
the Talbot pattern will be observable only within cer- 
tain boundaries; beyond these boundaries it will start 
to blur up, observing gradually the appearance of the 
Bragg diffraction channels. Instead of going to this limit 
by propagating quantum trajectories further and further 
away, let us rather consider the case where we increase 
the number of slits progressively until reaching a num- 
ber such that within a certain bound region we have the 
certainty of observing a Talbot pattern, but that further 
away we will only see Fraunhofer diffraction channels. 
This transition can be observed, from top to bottom, in 
Fig. [4] with both pt (left) and the associated quantum 
trajectories (right). 

As can be seen in Figs. IH^a) and Ufa'), the single slit 
case, the lack of neighboring trajectories allows the tra- 



jectories to spread out in all directions with no bound and 
therefore no (Talbot) pattern can be observed. Free prop- 
agation of a Gaussian wave packet just means, according 
to this trajectory picture, that one can observe free mo- 
tion (rectilinear and uniform) as in classical mechanics 
almost since the beginning of the propagation. In terms 
of the standard quantum mechanics, this means that pt 
spreads linearly in time, but keeping its volume constant. 
This constant volume only means that the number of tra- 
jectories conserves although they are far apart from one 
another. 

When two slits are considered, as in Figs. IHb) 
andUfb'), things change dramatically. From standard 
quantum-mechanical viewpoint we observe a channeling 
structure due to the interference of the two outgoing 
Gaussian wave packets. This is the Fraunhofer diffrac- 
tion pattern. From the quantum trajectory viewpoint, 
though trajectories leave along those different channels 
(with a very low density of trajectories in between), the 
fact that the perpendicular semiplane behind the slits 
has been divided into two identical halves, where trajec- 
tories do not cross the half one dominated by the opposite 
slit is remarkable. There is a very strong quantum pres- 
sure exerted by the trajectories arising from each half 
along the symmetry axis of the system. Moreover, it is 
also worth commenting that very close to the slits a cer- 
tain pattern, with two temporary maxima just behind 
the slits, is already present. In Figs. SJc) andlH^c'), the 
three slit case, a similar pattern to that described when 
discussing Fig.[3]^b) is observed. However, a Talbot pat- 
tern is still not formed because Fraunhofer diffraction 
channels emerge immediately. 

From the previous comments we can then establish 
that the Talbot effect can be seen as a "partition" of 
the space as the number of slits increases because of the 
strong effect (acting as an infinite barrier) of the quantum 
pressure. This is confirmed when we go to -/V = 10 [see 
Figs.Hd) andlHd')] andiV = 50 [see Figs. He) andHe')]. 
Indeed, the existence of the quantum pressure leads to a 
sort of quantum equilibrium state in which the unit struc- 
tures formed by the bunches of trajectories can coexist. 
Only when the quantum pressure starts decreasing, these 
units begin to blur up since the trajectories spread out 
the corresponding boundaries. Since this effect is similar 
to a dissipation, it could be called a quantum trajectory 
dissipation (though its nature is different to that of real 
dissipative phenomena). This "nonequilibrium" situa- 
tion remains until a new equilibrium is established: the 
Fraunhofer regime. One must realize that the stationar- 
ity of the Talbot regime is only typical of this near field 
phenomenon. In general, near field or Fresnel phenom- 
ena are not stationary. Nonetheless, as we will see in the 
next section, one can still speak about a certain class of 
stationarity (in the momentum space) within the Fres- 
nel diffraction regime different from that observed in a 
Talbot regime. 
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V. ATOM-SURFACE SCATTERING 

Talbot patterns appear whenever one deals with any 
kind of infinitely periodic structure in which interference 
can be observed, and not only with gratings consisting of 
slit arrays. This is the case, for instance, of atom-surface 
scattering where the surface plays the role of the periodic 
grating . ^'^'^^ In Fig.[5]the formation of the Talbot pattern 
is plotted when we have a beam of He atoms illuminat- 
ing ten unit cells of the Cu(llO) surface at perpendicu- 
lar incidence and energy Ez — 21 meV. This scattering 
system is described by a classical interaction potential 
V{x, z) = Vm{z) + Vc{x, z), where 



Vm{z) = D{1 - e-'^')^ - D 



(35) 



is a Morse potential, and 



Vcix, z)^De 



-2az 



0.03 cos 



2ttx 



0.0004 cos 



/ Anx 



(36) 



the coupling term between the two degrees of freedom. 
Here, D = 6.35 meV, a = 1.05 A"!, and d = 3.6 A (as 
the aperture of the slits used before) i^^'^'' The corruga- 
tion of the surface is very weak, as can be appreciated by 
the amplitudes of the cosine functions. 

In order to make the plot clear, we have neglected the 
incident part of the trajectories and considered as before 
the evolution of those corresponding to the centroid line 
of the incoming wave parallel to the surface. As can be 
seen in Fig.[5fa), the pattern is very similar to that dis- 
played by the ten-slit array plotted in Fig. [D(d'). This is 
because the corrugation of the surface is relatively weak. 
This can also be seen if we compare Figs.[5jb) and[3jb); 
though the latter has been obtained considering = 50, 
the structure is similar in both cases (we are inside the 
Talbot area). Note, however, that the Talbot structure 
is slightly distorted and repeats a bit further away than 
twice the Talbot distance considered before {2(P/X). The 
reason for this shift and distortion is the attractive part 
of the interaction potential V, which causes a certain 
acceleration in the motion of the particles. This is the 
so-called Beeby correction in atom-surface scatteringJ^ 
In other words, only when the potential is flat, we can 
assume that zt = (P/X. Otherwise, we should consider 
an effective Talbot distance, 



ZT 



X{x, z) 



ZT 



V{x,y) 



(37) 



where A — 2-Khl ^2m(Ez — y\ Taking into account 
that F < 0, it is clear that the square root factor in 
Eq. (|37p will be greater than 1, and therefore zt > zt- 
The correction factor in Eq. p7|) due to the presence 



of the well depth is y/l + D/E^, and therefore we ob- 
tain Zt I Zt = 1.14, which is basically the discrepancy 




FIG. 5: (a) Quantum trajectories corresponding to the He- 
Cu(llO) system at 21 meV and normal incidence. Only the 
emergent part has been plotted, (b) Enlargement of the tra- 
jectories shown in part (a). In both panels, the x distance is 
scaled in units of the Cu(llO) unit cell d, and z in units of 
twice the Talbot distance for an A^-slit grating {2zt ~ 2(f /X). 



observed in Fig. [IJb). Hence, within the context of 
atom-surface scattering, it would be appropriate to speak 
about the Talbot-Beeby effect, which gathers both the ef- 
fects caused by the periodicity (depending on Vc) and 
those arising from the attractive part of the interaction 
(depending on Vm)- 

It is clear that since there is no infinite beam of He 
atoms illuminating the Cu surface, the quantum pres- 
sure will decrease as the atoms get further away from the 
surface, and the Talbot pattern will disappear. However, 
there is something very interesting in this kind of sys- 
tems: one can extract already very important informa- 
tion about the diffraction peaks (experimentally detected 
at Fraunhofer distances) once the classical asymptotic re- 
gion is reached (V^ ~ 0). This is a nice manifestation of 
the effect mentioned above: the motion within the Fres- 
nel region is governed by the momenta that will give rise 
later on to the different Fraunhofer diffraction channels. 
This has been easily proven by Sanz et aL-21 by computing 
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FIG. 6: DifFraction of a particle with a mass 500 times heavier 
than that of a He atom by (a) ten slits and (b) ten unit cells 
of the Cu(llO) surface. In both panels, the x distance is 
scaled in units of the Cu(llO) unit cell d, and z in units of 
twice the Talbot distance for an A^-slit grating {2zt ~ 2d^/A). 
Moreover, the perpendicular incidence energy is Ez — 21 meV 
in both cases. 

the S'-matrix elements in the classical asymptotic region 
and comparing them with the intensity calculated from 
quantum trajectories collected in the asymptotic region. 
A similar calculation for slit arrays by the same authors 
can be seen in Ref. [l^. The Fraunhofer regime is reached 
very far from the region where the interaction potential 
V is negligible and increases its distance as the num- 
ber of unit cells illuminated by the initial atomic beam 
is increased. Thus, for example, for ten unit cells the 
Fraunhofer regime is reached around 1000 A. 



VI. CORRESPONDENCE PRINCIPLE AND 
TALBOT EFFECT 

As heavier impinging particles are considered, i.e., 
when we approach the classical limit, the behaviors ob- 
served in previous sections undergo dramatic changes. 



approaching or resembling those expected classically. If 
instead of He atoms we consider a fictitious particle with 
a mass 500 times that of a He atom, it is clear that 
quantum effects should disappear or, at least, decrease. 
This is what one could think, mistakenly, when observ- 
ing Fig.[S]Ja), where we have represented the diffraction of 
such particles by a grating consisting of ten slits. As can 
be seen, the trajectories are basically straight lines, what 
induces to think that no interference effects are present 
(neither Talbot nor Fraunhofer ones), but particles move 
tracing a simple rectilinear, uniform motion. Obviously, 
this is misleading; both the x and z directions are given 
in terms of the old d and zt- However, if we replot this 
figure taking into account that now zt has increased by 
a factor y^m/mB_c (where m is the mass of the fictitious 
particle), we obtain again a Talbot pattern. And far 
away we would find the same Fraunhofer pattern as be- 
fore. What has happened is that particles move now more 
slowly, and therefore the spreading of the corresponding 
bunches of trajectories will also be slower. 

The previous result shows that quantum particles re- 
main quantum even in the so-called classical limit;^^ one 
only needs to be patient and wait enough time in order to 
observe again the quantum phenomena (of course, unless 
the complexity of the system is so enormous that any 
quantum effect is imperceptible experimentally). How- 
ever, what happens if instead of ten slits we illuminate 
ten cells of the Cu(llO) surface with a beam of our fic- 
titious particle? The answer appears in Fig. [Bl^b): now 
the topology displayed by the quantum trajectories tries 
to resemble that of the classical one, but the noncrossing 
principle holds. That is, since classical trajectories give 
rise to the appearance of two caustics (direction of maxi- 
mum refiected intensity) at the so-called rainbow angles, 
the quantum trajectories will try to describe a similar 
structure, with the particularity that they cannot cross, 
and therefore, after tracing an almost straight line, they 
will be bounced backwards in relatively sharp angle, as 
seen in the lowest part of Fig. [Hb). Moreover, notice 
that, because of this motion, the laminarity of the fiow 
described by the trajectories is lost; now the topology of 
the trajectories is more irregular, this leading to some 
crossings (though they occur at different times). 



VII. CONCLUSIONS 

To conclude, there are three important points in this 
work that are worth stressing. First, regardless of the 
gradual spreading of the partial wave function associated 
with each unit cell of a periodic grating, the correspond- 
ing initial swarms of quantum trajectories keep moving in 
a bound space region. To some extent, this is analogous 
to considering the trajectories moving inside a multimode 
cavity.— Because of this effect, within the Talbot in- 
terferometry context, the superposition principle and the 
Born-von Karman boundary conditions in Bohmian me- 
chanics have to be understood in a very different way as 
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they are in standard quantum mechanics. It is not simply 
the overlapping of a number of wave functions but alloca- 
tion of identical copies of the same trajectory behavioral 
pattern. Second, this behavior is somehow similar to the 
classical behavior observed when studying periodic struc- 
tures: all the information about the periodic structure 
can be obtained by simply studying the classical effects 
provoked by one of the periods although nonlocal effects 
are always presentJ^ And third, the Talbot-Beeby effect 
is proposed to understand quantum trajectory carpets in 



the presence of attractive interaction potentials. 
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